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Abstract 

We upgrade the study of the physical reach of a Neutrino Factory in the Four 
Family Neutrino Mixing scenario taking into account the latest LSND results 
that points out how the 3+1 scheme cannot be completely ruled out within the 
present experimental data (although the 2+2 scheme is still the preferred choice 
when four neutrinos are considered) . A detailed comparison of the physical reach 
of the i^-factory in the two schemes is given, with similar results for the sensitivity 
to the mixing angles. Huge CP-violating effects can be observed in both schemes 
with a near, O(10) Km, detector of O(10) Kton size in the Vn — * v T channel. A 
smaller detector of 1 Kton size can still observe very large effects in this channel. 
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1 Introduction 



Indications in favour of neutrino oscillations have been obtained both in solar neutrino 
0, |, [| £|, [| and atmospheric neutrino || [7|, || || 0] experiments. The latest atmo- 
spheric neutrino data imply Am^ m ~ (1.6 — 4) x 10 -3 eV 2 [fTT]] , whereas the solar 
neutrino data prefer Am 2 nn ~ 10~ 10 or 10 -7 — 10 -4 eV 2 , depending on the particular 
solution for the solar neutrino deficit. The LSND data [1^, [13]], on the other hand, 
would indicate a — > v e oscillation with a third, very distinct, neutrino mass differ- 
ence: Am 2 LSND ~ 0.3 — 6 eV 2 . The LSND evidence in favour of neutrino oscillation 
has not been confirmed by other experiments so far. However, the MiniBooNE exper- 
iment Jl4| will be able to confirm it or not in the near future. If the LSND results are 
confirmed we would face three independent evidence for neutrino oscillations character- 
ized by squared mass differences quite well separated. To explain the whole ensemble 
of data at least four different light neutrino species are needed. The new light neutrino 
is denoted as sterile [15], since it must be an electroweak singlet to comply with the 
strong bounds on the Z° invisible decay width |nj. We stress that three massive light 
neutrinos can not explain all the present experimental results, as it has been shown 



with detailed calculations in [|T?J . 

There are two, very different, classes of four neutrino spectra: three almost degen- 
erate neutrinos and an isolated fourth one, or two pairs of almost degenerate neutrinos 
divided by the large LSND mass gap. The two classes of mass spectra are usually 



called the 3+1 and 2+2 schemes |L8|] , respectively. All the present experimental evi- 
dence for neutrino oscillations have been combined in the literature in order to identify 
which of the two classes of mass spectra better adapts to the data. The experimental 
results were strongly in favour of the 2+2 scheme [[19] until the latest LSND results 
have been presented in June 2000, p(J (see [13]). The new analysis of the experimental 



data results in a shift of the allowed region towards smaller values of the mixing angle, 
sin 2 (29 lsnd), reconciling the 3+1 scheme with exclusion bounds coming from CDHS 
0, CCFR (22| and Bugey |Z3J. Although the 2+2 scheme is still favoured^, the 3+1 
scheme is at present marginally compatible with the data, [^, [27], |28[] . However, the 
2+2 and the 3+1 scheme face the upcoming experiments on totally different footing: if 
MiniBooNE disconfirms LSND, the 2+2 scheme is falsified. On the contrary, it is not 
possible to falsify the 3+1 scheme: we can always consider an extension of the Standard 
Model with three light neutrinos and a fourth sterile one, separated by some squared 
mass difference, Am^ 2 3)-4- The implication of a negative result of MiniBooNE is just 

The specific form of the neutrino mass spectrum appears, therefore, one of the 
(many) open questions related to the lepton sector of the Standard Model. 

Four neutrino oscillations imply a Maki-Nakagawa-Sakata (MNS) 4x4 mixing ma- 



3 A novel Bayesian analysis of the exclusion bounds, in the spirit of |l9]] has been presented in 
claiming that the 3+1 scheme is allowed at the 99 % CL only, but not at the 95 % CL. 



2 



trix, with six rotation angles % and three phases 8{ (for Majorana neutrinos, three 
additional phases are allowed, but they are not testable in oscillation experiments, and 
therefore will not be considered here). A Neutrino Factory ||29| , |30H seems to be the best 
option to explore this huge parameter space. The /i ± -decay into the straight section 
of a muon storage ring should produce a very intense and pure neutrino beam. The 
rich flavour content (50 % of v^ipy) and 50 % of z/ e (z/ e ) are simultaneously produced), 
finally, makes the Neutrino Factory well suited for precision studies of the MNS mix- 



ing matrix, hopefully including the discovery of leptonic CP- violation [[31], R^]. The 
following scheme reminds that at the Neutrino Factory \x and r appearance channels 
can also be used, in combination with the [i and e disappearance experiments. 



fi -»• e Vp v e -> (i ,e T (1) 
i 

i 

In |33| , [341 , |35| the "wrong-sign muon" channel (/x + appearance in a fi~ beam) has been 
shown to be extremely useful to explore the parameter space of three-family neutrino 
mixing, with particular interest in the measure of the (single) CP-violating phase, thus 
deserving the nickname of "golden measurement" at a Neutrino Factory. This has to 
be compared with a conventional beam experiment (using muon neutrinos from pion 
decay), such as K2K or the approved FermiLab to Soudan long baseline experiment. 
In these experiments, mainly the //-disappearance channel is exploited. 



In [gg, [37] ^ was shown that a Neutrino Factory with 10 — 50 GeV muons can attain 
sin 2 (0y) as low as 10~ 5 - 10~ 3 for Am 2 LSND e [10 _1 , 10 1 ] eV 2 . Moreover, it was found 
that sizeable CP-violating effects can be observed in the — *> v T channel with a 1 Kton 
detector located at O(10) Km. This analysis has been performed in the 2+2 scheme, 
the only allowed at that moment. The first motivation for this paper is, therefore, 
the comparison of what has been found for the 2+2 scheme with the same kind of 
analysis in the, by now marginally allowed, 3+1 scheme. A better understanding of 
the oscillation probability structure is a natural by-product of this analysis, both for the 
CP-conserving and the CP-violating part. In particular, we found a simple argument 
that shows how the — > v T channel is the best suited one for CP-violation experiments 
in four neutrino mixing, to be compared with the three-family mixing where v e — > 
happens to be optimal. Finally, the same argument justifies the loss in sensitivity to 
small mixing angles in the 3+1 scheme with respect to the 2+2 scheme. 

The considered set-up is, as in |36], [37I , a Neutrino Factory with 2 x 10 20 yU + and \i~ 
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decaying in the straight section of a 10 — 50 GeV muon storage ring per year, and five 
years of data taking. Muon energies in this range are at present under discussion. The 
higher energy range allows a good background rejection ||381 ; moreover, the integrated 



flux times the cross-section increase with E^. A high-energy Neutrino Factory seems 
therefore the best option, with the energy mainly limited by cost considerations. How- 
ever, although the total number of charged leptons into the detector increases with the 
parent muon energy, the flux of low energy neutrinos decreases. If low energy neutri- 
nos are needed (for example, to study CP-violating observables strongly reducing the 
matter effects |39|, HUH), this reduction in the flux should be taken into account. 



The mixing angles that relate neutrino mass eigenstates with an LSND mass dif- 
ference can be studied in short baseline experiments, L ~ 1 Km. A small detector 
with 0(1) ton mass is, therefore, well suited to study the whole gap-crossing param- 
eter space (due to the large neutrino flux that illuminates the detector). To take full 
advantage of the rich flavour content of the beam, this detector should be equipped 
with r-tracking and (/i, r) charge identification. If CP-violating observables are con- 
sidered, the best option is a larger, O(10) Kton detector located at O(10) Km down 
the source. This set-up is equally powerful both for the 2+2 and the 3+1 schemes, and 
should be compared with the typical set-up needed when three-family neutrino mixing 
is considered. 

We also try to answer to the following question: is it possible to explore the whole 
parameter space with a different detector, with no r-tracking, but taking full advantage 
of the energy dependence of the transition probabilities, in the spirit of We focus 
on the v e — > Un channel in the 3+1 scheme, with a realistic 10 Kton magnetized iron 



detector of the type presented in |38| located at L = 40 Km down the source. The 
Neutrino Factory is run with 2 x 10 20 useful muons per year for 5 operational years 
for both muon polarities, at = 50 GeV , with a detector energy resolution of 
AE U = 10 GeV . A detailed estimate of the backgrounds and detection efficiencies 
of the considered detector has been presented in |33|. The energy dependence of the 



oscillation probabilities could in principle help in the measurement of two gap-crossing 
angles, or one angle and a CP-violating phase, at a time. In the latter case, we find 
results similar to those in [Q: we can easily reconstruct the phase and the angle at 



the same time, with an error of tens of degrees on Si and of tenths of degree on the 
angle. However, it seems extremely difficult to measure two gap-crossing angles at a 
time in the v e — > channel. Our conclusion is that to fully explore the parameter 
space of the four-family model a detector with r-tracking is needed. 

The paper is organized as follows: in Sect. |^ we introduce our parametrization of 
the MNS mixing matrix both for the 2+2 and 3+1 schemes; in Sect. ^.1| the present 
bounds on the mixing angles coming from existing experiments are given; in Sect. |2]^ 
we describe the Neutrino Factory and detector setup; in Sect. ^ we present our results 
for the sensitivity of the Neutrino Factory to the mixing angles (in the case of no CP 
violation), comparing the 2+2 and 3+1 schemes; in Sect. |] we extend our analysis to 
the CP-violating observables; in Sect. |5|we study the possibility of measuring two gap- 
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crossing angles or one angle and a phase at a time exploiting the energy dependence 
of the transition probabilities; in Sect. ^|we eventually draw our conclusions. 



2 The Four Neutrino mixing matrix 

When four neutrinos are considered, two very different classes of mass spectrum are 
possible: three almost degenerate (mainly active) neutrinos, accounting for the solar 
and atmospheric oscillations, separated from the fourth (mainly sterile) one by the large 
LSND mass difference, Am 2 LSND ; or, two almost degenerate neutrino pairs, accounting 
respectively for the solar and atmospheric oscillations, separated by the LSND mass 
gap. The two mass spectrum classes are depicted in Fig. [TJ We refer to these possibil- 
ities as 3+1 and 2+2 scenarios. There are four 3+1 and two 2+2 scenarios depending 
on the specific ordering of the mass differences. Notice that the intriguing hierarchical 
and inverted hierarchical mass spectrum are 3+1 scenarios. 



Class I 



A 

Airu 



Class II 



B 



Arn 



Am, 



Am 



Am , 



Am , 



Am' 



Am, 



Figure 1: Different types of four family neutrino mass spectrum: 3+1 scenarios (left); 
2+2 scenarios (right). 



It has been shown in |L9| that the combined analysis of solar, atmospheric and LSND 
data disfavours the 3+1 scheme. For this reason, the 2+2 scheme has been carefully 
studied in the recent literature (see, for example, |^T], [42] and references therein). 
Consider for definiteness the rightmost scenario, namely the lower pair accounting 
for the solar neutrino deficit and the higher pair for the atmospheric one (the other 
possibility directly follows by changing the sign of the Am 2 LSND ). The is therefore 
in the heavier pair and the v e in the lighter one. Is the sterile neutrino, z/ s , responsible 
for the observed atmospheric oscillations or for the solar neutrino deficit? The latest 
SuperKamiokande results for the atmospheric neutrinos^ disfavours total conversion 

4 In particular, the zenith angle distribution of the upward-going muons p3[ | , the partially contained 
multi-ring events and the neutral current data sample. 
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of Ufj, into the sterile neutrino at 99 % CL ||44|| . Moreover, the conversion of u e solar 



neutrinos into active neutrinos gives better global fits of the experimental data, with 
respect to active-into-sterile conversion ||45|| . Although partial conversion of — > v s 



or v e — > v s is not excluded (with a sterile component in the atmospheric oscillation as 
large as 50 %, f46fl), the present solar and atmospheric data suggest active-to-active 
oscillations. The active-to-sterile oscillation should therefore be responsible only for 
the LSND results. This scenario appears quite unnatural in the framework of the 2+2 
scheme. 



The latest analysis of the LSND data f20fl , however, shows a shift of the allowed 
region for the LSND two-family-equivalent mixing angle, sin 2 (29 lsnd), towards smaller 
values. This reconciles the 3+1 scheme with the exclusion bounds coming from CDHS 



2jJ , CCFR [|22j and Bugey In the 3+1 scheme, the three almost degenerate 

neutrinos are mainly active and the separated fourth is mainly sterile; the gap-crossing 
mixing angles are generally small. In this scenario, the interpretation of solar and 
atmospheric oscillations as active-to-active and LSND as active-to-sterile naturally 
arises: this scheme is a deformation of the Standard Model, slowly decoupling as the 
gap-crossing mixing angles become smaller and smaller (and thus, the sterile becomes 
irrelevant). 

Given n light neutrino species, the most general mixing matrix is an n x n unitary 



matrix ||47|| , Umns- For n = 4, the MNS matrix contains six independent rotation 
angles 6ij and three (if neutrinos are Dirac fermions) or six (if neutrinos are Majorana 
fermions) phases <5j. However, oscillation experiments are only sensitive to the first 
three phases, the effect of the Majorana phases being suppressed by factors of m v /E v . 
The Majorana or Dirac nature of neutrinos can thus be tested only in AL = 2 tran- 
sitions such as neutrino-less double /5-decay ||4~8| . In the following analysis, with no 
loss in generality, we will restrict ourselves to Dirac-type neutrinos only. We consider 
a hierarchical 3+1 spectrum and a class II-B 2+2 spectrum, for definiteness. 

This large parameter space (6 angles and 3 phases, to be compared with the stan- 
dard three-family mixing case of 3 angles and 1 phase) is actually reduced to a smaller 
subspace whenever some of the mass differences become negligible. Consider the mea- 
sured hierarchy in the mass differences, 

Am 2 sol < Ami < ^rn 2 L SND , (2) 

and define 

At short distance, L = 0(1) Km, for neutrinos up to 0(10) GeV, 

A so ; ) A atm <C 1 , 
A L snd = 0(1) . (4) 

Therefore, it is natural at short distances to neglect the solar and atmospheric mass 
difference and to work in a reduced parameter space. This approximation is called 



6 



"one-mass dominance" [ 4"5| . In the 2+2 scheme, neglecting the smaller mass differ- 
ences implies that rotations in the (1 — 2) and (3 — 4) planes are irrelevant. Thus, it is 
not possible to measure the rotation angles in these planes in oscillation experiments. 
Two CP- violating phases also become irrelevant, and therefore the reduced parameter 
space in the 2+2 scheme contains 4 rotation angles and 1 phase. In the 3+1 scheme, ne- 
glecting the solar and atmospheric mass differences implies that rotations in the whole 
three-dimensional subspace (1 — 2 — 3) are irrelevant for oscillation experiments, and 
the physical parameter space contains just three rotation angles and no phases. When 
considering CP-violating phenomena, however, at least two mass differences should 
be taken into account: in this case we neglect the solar mass difference and consider 
the atmospheric mass difference a perturbation. This is called "two-mass dominance" 
approximation. In this approximation, regardless of the scheme, the parameter space 
contains 5 angles and 2 phases. The number of independent parameters of the MNS 
mixing matrix in four neutrino models is summarized in Tab. |l[ 





Angles 


Dirac CP-phases 


Majorana CP-phases 


Major ana z/s 


6 


3 


3 


Dirac z/s 


6 


3 





Dirac i/'s 
Amj 2 = 


5 


2 





2+2 
Dirac z/s 
Amf 2 = ATO34 = 


4 


1 





3+1 
Dirac z/s 
Amf 2 = Am^ 4 = 


3 









Table 1: Parameter space in four neutrino models: for Dirac neutrinos we consider the 
general case (three non-zero mass differences) and the one- and two-mass dominance 
approximations; for Majorana neutrinos we consider the general case only. 



A generic rotation in a four dimensional space can be obtained by performing six 
different rotations in the ii — j) plane, resulting in plenty of different parametriza- 
tions of the mixing matrix (and still not taking into account the three CP-violating 



7 



phases). However, in |36| , |37|| was shown how the one-mass dominance and two- mass 
dominance approximations can be implemented in a trasparent way (in the sense that 
only the physical parameters appear in the CP-conserving and CP- violating oscillation 
probabilities). A convenient parametrization of the mixing matrix is that in which the 
rotation matrices corresponding to the most degenerate pairs of eigenstates are located 
at the extreme right. If the eigenstates i and j are degenerate and the matrix Uij is the 
rightmost one, the corresponding angle Oij automatically disappears from the oscilla- 
tion probabilities, and the parameter space gets reduced to the truly observable angles 
and phases. If a different ordering of the rotation matrices is taken, no angle disap- 
pears from the oscillation formulae, and a parameter redefinition would be necessary 
to reduce the parameter space to the observable sector. 



In the 2+2 scheme, the following parametrization was adopted in |56| implementing 
the previous argument: 

Umns = U u (9u) U 13 (9 l3 ) U 24 (6 24 ) U 23 (9 23 , 5 3 ) U 34 (6 34 , 8 2 ) U l2 {6 12 , 5 X ). (5) 

In the one-mass dominance approximation, the unphysical angles and phases (6i 2 ,5i) 
and (6*34, 5 2 ) automatically decouple. The oscillation probabilities in the appearance 
channels areQ: 

Plf{v e -> up) = 4^444 sin 2 (^^] , (6) 



Pl¥(y e ^v T ) = 444 



AE 

(2 2 2,2 2 2 \ 
\ S 13 S 14 S 23 ' C 14 C 23 S 24J 



2 f Amn 3 L\ 

2C14S14C23S23S13S24 cos d 3 \ sin — 7^— , (7) 



AE 



p2+2( \ a 2 2 1 I 'I 'I , 'I 'I 'I \ 

^CP \ V H ~^ V t) — 4C 23 C 13 l S 13 S 14 C 23 + C 14 S 23 S 24J 



,2 2 2 1 2 2 „2 



+ 2ci 4 Si 4 C23S23Sl3S24COS()3 Sin 1 



AE 



in the disappearance channels are: 



Pcp>,^^) = 1 - 444(4, + 44) sin 2 (^f^ ), (9) 
- f e ) = i - 444(4 + 44) sin 2 ( ■ (10) 



CP 



\ AE 



Notice that the physical phase £3 appears in the CP-conserving transition probabilities 
in a pure cosine dependence. No CP-odd observable can be built out of the oscillation 
probabilities in this approximation in spite of the existence of a physical phase in the 
mixing matrix. 



5 In what follows, we separate the CP-even terms from the CP-odd ones: 
P(y a -> vp) = P C p iy a ™> vp) + Pofr (u a -> vp). 
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In the two-mass dominance approximation, new CP-violating terms arise. Expand- 
ing the probabilities at first order in A atm , we getf]: 



Pop iye Vy) = -8c 13 c 23 C24C34S23S24.S34 sm(S 2 + S 3 ) 
I Am| 4 L\ . 2 / Am 23 L 



sin' — ^ , 11 



AE V ) \ 4E, 



PQ?{y e -> V T ) = 8C 23 C 24 {c23C 34 S23S24S34(Ci4 - S^S^) sin(5 2 + h 



+ C14C34S13S14S34 



'24 



s 23 ) sin <5 2 — s 23 s 24 sin(^2 + 28. 



+ Ci4C24Si3Si4S23S24(c 34 — S 34 ) Sm5 3 | 



^ ^t) = -8ci 3 c 23 Ci 4 c 24 c 3 4S 3 4 [c 23 si 3 si4 sin5 2 + C14S23S24 sin(<5 2 + 5 3 )} 

x (^) (^) ■ <13) 

Two distinct phases, 52 and £3, appear in these formulae in a characteristic sine de- 
pendence which is the trademark of CP-violating observables. CP-violating effects can 
only be measured in appearance channels, whereas the disappearance channels u e — > v e 
and Vy — >■ z/^ are only sensitive to the CP-even parameters, 

P&p{v e -> f e ) = P&p{Vy -> ^) = P&p{v T -»• «r) = ■ (14) 

In the 3+1 scheme, the following parametrization shares the same virtues of eq. fl5|): 

£W = MM 024(024) 034(^34) 023^23 , 01 3 (013 , <* 2 ) 01 2 (012 , fc). (15) 

This parametrization has the additional advantage that the three-family model mixing 
matrix in its standard form can be immediately recovered when 9^ = 0. For small 
gap-crossing angles 8^, we expect slight modification with respect to the three-family 
model. 

In the one-mass dominance approximation, the unphysical angles and phases (6i 2 , Si), 
(013) S 2 ) and (023,^) automatically decouple. The oscillation probabilities in the ap- 
pearance channels are: 

Pl\\ve^v») = 4^444 sin 2 f^f^J , (16) 

P&\v, - u T ) = 444^4 sin 2 (^f^) , (17) 

Pcpiyn^Vr) = 4c 2 4 44 sin 2 ( Am l 4L ) ; (18) 



4E 



3 In |3q, p7[ some misprints were present in the v e — > v T formula that have been corrected here. 
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in the disappearance channels are: 

Pht\v^v») = 1 -444(4 + 44) sin 2 {^f) , (19) 

Pl\\v e ^v e ) = i -4444(1 -444) sin 2 [^^] ■ (20) 

As already stressed, angles and phases in the three-dimensional physically irrelevant 
subspace are not present in these formulae. 

Finally, in the two-mass dominance approximation we get, expanding at first order 
in A atm : 

P&piVe^V^) = 8 4 C 23C24Sl4S 2 4 {-Cl3Si 4 S23S34Sin5 2 

+ C14S13 [C13C23S24S34 sin S 3 + C24S23 sin(5 2 - 5 3 )]} 
f Aml 3 L\ . 2 { Am 2 u L\ 

x hiH sm Mr)' (21) 

P&p{v e -> v T ) = 8 4 c i3C23C 2 4Si 4 s 3 4 [s14.s23.s24 sin 5 2 - C14C23S13 sin 83] 

x hf-J 81 " 2 hf-J • (22) 

PQptyp -> V T ) = -8 4 C 23C 2 4S2 3 S24S 3 4 SU1 8 2 

x (-if-J^ i-if-J- (23) 

These formulae, both in the 2+2 and the 3+1 scheme, will be used in Sect. § 
and Sect. |] to explore the parameter space of the four-family model at the Neutrino 
Factory. 



2.1 Experimental bounds on the gap-crossing angles 

We recall here the bounds on the rotation angles and mass differences coming from the 
existing experiments. The Bugey and CHOOZ experiments P3L BUf give strong upper 



limit to the v e — > v e disappearance two-family equivalent mixing angle. In two families, 
P CP (ye - ve) = 1 - sm\26) exp sin 2 ( ^<sndL ^ (m) 

with sm 2 (29) exp < 0.2 in the LSND-allowed region. The positive result from LSND 
gives a lower limit on the u e — > two-family-equivalent mixing angle, 

Pcp(^ - v e ) = sm 2 (26) LSND sin 2 (25) 
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with 10 3 < sin 2 (29)LSND < 1- These bounds, jointly with the negative results from 
Karmen2 



must be 



JlJ] and previous experiments such as CDHS and CCFR |2l], [ 

interpreted in the 2+2 and 3+1 scheme, extracting informations on the gap-crossing 
angles and mass differences. 

• The 2+2 scheme 

In the 2+2 scheme, that is still favoured by the data, the bound on v e disappear- 
ance translates into an upper limit on the combination 



4, sin 2 (2# 24 ) +4 sin 2 (29 23 
whereas the bound on v e appearance implies 

10~ 3 < c 2 3 c4sin 2 (2fl ; 



< 0.2 



723 J 



< 10" 



(26) 



(27) 



These bounds suggest the conservative (or even "pessimistic" ) hypothesis adopted 
in P5| : to consider the four gap-crossing angles 6*13, 6*14, #23 and #24 to be equally 
small in the mass difference region Am 2 LSND e [10 _1 , 10 1 ] eV 2 . We follow here 
the same hypothesis: all the gap-crossing angles are small (i.e. less than 10°), with 
the possible exception of one angle that we leave to vary in some interval. The 
remaining angles 612 and #34 are directly the solar and atmospheric mixing angles 
in the two-family parametrization, respectively. The typical flavour content of 
the mass eigenstates in the 2+2 scheme is presented in Fig. |2| (left). 

The 3+1 scheme 

This scheme is only marginally allowed (a recent study shows that it is 
compatible with the experimental data at the 99 % CL only). However, it is a 
natural extension of the three-family model. There are four very small allowed 



region in the two-family equivalent parameter space |25| : 



1. 


Am 2 4 


~ 0.3 


eV 2 


sin 2 (26») L 5Ar D 


~ 2 x 10 


2. 


Amf 4 


~ 0.9 


eV 2 


sin 2 (26») L5A r D 


~ 2 x 10 


3. 


Amf 4 


~ 1.7 


eV 2 


sin 2 (26») L 5jv D 


~ 1 x 10 


4. 


Amg 4 


~ 6.0 


eV 2 


sin 2 (29) LS nd 


~ 2 x 10 



We restrict ourselves to case 2, for simplicity. In this case, we get in our 



parametrization for the v e appearance mixing parameter 

4 A s\ 4 sin 2 (2^4) ^ 2 x 10" 3 . 



f28l 



This bound is consistent with the conservative hypothesis of equally small gap- 
crossing angles that will be followed in the rest of the paper. In the 3+1 
scheme the remaining angles, 6*12, 6*23 and #13 can be obtained by the combined 
analysis of solar and atmospheric data in the three-family parametrization. The 
typical flavour content of the mass eigenstates in the 3+1 scheme is presented in 
Fig. I (right). 
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Figure 2: The flavour content in the mass eigenstates with a representative choice for 
the mixing angles: in the 2+2 scheme with 612 = 45°, #34 = 45°, #13 = #14 = #23 = #24 = 
5° (left); in the 3+1 scheme with 12 = 45°, 13 = 13°,0 23 = 45°, i4 = 6 24 = 6 U = 5° 
(right). The different flavours are, from lightest to darkest: v s ; v,; v e and v T . 

2.2 Experimental Setup: The Neutrino Factory and the De- 
tector 



In the muon rest frame, the distribution of muon antineutrinos (neutrinos) and electron 



neutrinos (antineutrinos) in the decay ji 



e ± + u e ( 



is given by: 



d 2 N 
dxdVL 



^[/o(^)T^A(x) costf] 



(29) 



where E v denotes the neutrino energy, x = 2E v jm iX and V, is the average muon 
polarization along the beam directions. $ is the angle between the neutrino momentum 
vector and the muon spin direction and m M is the muon mass. The positron (electron) 
neutrino flux is identical in form to that for muon neutrinos (antineutrinos), when the 
electron mass is neglected. The functions fo and fi are given in Table |2|, |52] . In the 
laboratory frame, the neutrino fluxes, boosted along the muon momentum vector, are 
given by: 



d 2 N D 



dydVt 



d 2 N„ 



dydVt 



4n, 



Ey (1-/3 cos <p) { 



3m, 



(3 cos if) 



TV, 
24n M 



m. 



4E£y (1-/3 cos p)]} 



(3 cos ip) I 



m, 



2E\y{\ -/3 cos tp) 



T^[mJ-2^(l-/3cos^)]} 



(30) 





fo(x) 


fi(x) 




2x 2 (3-2x) 


2x 2 {l- 2x) 




12x 2 (l - x) 


I2x 2 (l - x) 



Table 2: Flux functions. 
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Here, (3 — J 1 — m 2 /i? 2 , E^ is the parent muon energy, y = E u /E^, n M is the number 
of useful muons per year obtained from the storage ring and L is the distance to the 
detector, tp is the angle between the beam axis and the direction pointing towards the 
detector. We shall consider in what follows as a "reference set-up" a neutrino beam 
resulting from the decay of n M = 2 x 10 20 unpolarized positive and/or negative muons 
in one of the straight sections of a muon storage ring (i.e. we do not consider two 
baselines operating at the same time) per year. The collected muons have energy E^ 
in the range 10 — 50 GeV. This energy range is under discussion as a convenient goal 
(a definite answer on which is the optimal energy to run the Neutrino Factory is still 
missing). The angular divergence 5<p is taken to be constant, 5p ~ 0.1 mr. 

The charged current neutrino and antineutrino interaction rates can be computed 
using the approximate expressions for the neutrino-nucleon cross sections on an isoscalar 
targetf], 

E E 
a uN ~ 0.67 x 10~ 42 x x m 2 a pN ~ 0.34 x 10" 42 x x m 2 . (31) 

To explore the whole CP-conserving parameter space we need a detector with r 
tracking and (/i, r) charge identification capability. As the dominant signals are ex- 
pected to peak at L/E v ~ 1/ Am 2 LSND , most of the parameter space can be explored 
in short baseline experiments (SLB), with L ~ 1 Km. At such a short distance from 
the source the neutrino flux is so intense that a small detector is well suited to study 
CP-conserving transitions. In what follows we consider an hypothetical 1 ton detector, 
with no detailed calculation of background and efficiencies as a function of the neutrino 
energy. We consider a constant background B at the level of 10 -5 of the expected num- 
ber of charged current events, Ncc, and a constant reconstruction efficiency e M = 0.5 
for and e r = 0.35 for r ± . The number of expected charged leptons in absence of 
oscillation is N^- = 9.3 x 10 8 and N e + = 4.0 x 10 8 for a fi~ beam ( = 4.7 x 10 8 and 
N e - = 7.9 x 10 8 for a fi + beam). We also applied a conservative cut on the neutrino 
energy: neutrinos with E v < 5 GeV have not been included in our results. 

To extend our analysis to the CP- violating parameter space, a larger detector must 
be considered: we choose an hypothetical 10 Kton detector, located a bit farther from 
the neutrino source, at L = O(10 — 100) Km. Also in this case fi and r charge 
identification capability is needed, being the — * v T transitions the optimal channel 
to observe CP- violation in a four-family model (as will be explained in the following). 
The same background and reconstruction efficiencies as for the CP-conserving sector 
are included, and again a fiducial cut on neutrinos with E u < 5 GeV is applied. 



7 For the i/ r -nucleon interaction, we used the exact expression for the cross section taking into 
account the r-mass. 
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3 Sensitivity reach of the Neutrino Factory 



We concentrate now on the sensitivity to the different gap-crossing angles that appear in 
the oscillation probabilities when only the LSND mass difference is taken into account, 
namely eqs. (|6|-|T0|) for the 2+2 scheme and eqs. ([T6| - p0] ) for the 3+1 scheme^. 



We define the sensitivity in the appearance channel as follows: the number of total 
expected events for a given flavour p a is 

N tot = N* ± AN* + Np (32) 

where 

= N a -B, (33) 

N p = N a < P{y a ^Up)>, (34) 

with N a the number of expected events in the absence of oscillation, B the fractional 
background (we consider B = 10~ 5 ) and < P[y a — > i/p) > the transition probability 
averaged over the v a flux and the CC interaction cross-section. Fluctuations over the 



background are taken to be gaussian, AN® = yiVj*. The excluded zone at 90 % CL 
(following f53| ) if no event is observed is the region to the right of 



Np = 1.65 AN" . (35) 

The sensitivity in the disappearance channel is defined as follows: the number of total 
expected events for a given flavour u a is 

N tot = N a -(l-B)± A[N a ■ (1 - B)] - N p (36) 

where 

Np=J2 N <*< p ("a -> Up) > (37) 

summing over all flavours distinct from v a . In this case we compare Np with the 
gaussian fluctuation over N a x (1 — B) (we notice that a background B at the level 
of 10~ 5 plays a marginal role, with respect to the appearance case). Again, following 
if no event is observed the region to the right of 



Np = 1.65 A[N a ■ (1 - B)} (3S 

is excluded at 90 % CL. 



8 Although we refer to the one- or two-mass dominance approximation (in the section devoted 
to CP-violating observables) formulae, all the numerical results have been obtained with the exact 
expressions for the transition probabilities. 
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3.1 Sensitivity in the 2+2 scheme 



We recall here the results of |36|, [37], albeit rederived with slightly different input 
parameters. In the one-mass dominance approximation, the CP-conserving parameter 
space consists of four rotation angles (#i 3 , 6*14, 6*23 and 6*24) and one phase, 5 3 . In what 
follows we set £3 = 0. The useful channels to measure or put severe upper limits on 
the gap-crossing angles at the Neutrino Factory are the following (for a fi~ decay): 



Ve - 








' appearance) 


Vp ■ 




-> \T 




" disappearance) 


Ve ■ 


-> v T - 






appearance) 




-> v T - 






appearance). 



(39) 

In order to present the sensitivity to a specific sin 2 9, we adopt the following ap- 
proach: we vary sin 2 9 between 10~ 7 and 1; the remaining three angles are considered 
to be already known: two of them are fixed to a small value, 9ij = 2°, and the third 
one is varied from 1° to 60°. The remaining parameters (those measured in solar and 
atmospheric experiments) are taken as follows: 

012 = 45°, #34 = 45°; 
Am 2 2 = l(r 4 eV 2 , Am 2 M = 3.5 x l(T 3 eV 2 . 

The large mass difference Am| 3 is varied from 10~ 3 to 10 2 eV 2 . At L — 1 Km matter 
effects are not relevant, since such a baseline is short enough to be completely above 
ground. We consider 2 x 10 20 useful muons per year and 5 years of data taking, with 
En = 20 GeV . For simplicity, the Neutrino Factory is supposed to be working with 
negative muons only. 

• Sensitivity to sin 2 #23 : ^ + appearance 

The /i + appearance channel (the so-called "wrong-sign" muons) is particularly sen- 
sitive to #23- Fig. |3] shows the 90 % CL exclusion curve in the sin 2 6*23 / Am^ plane 
for different values of 613. In the LSND-allowed region (Ami, in the lO" 1 - 10 1 eV 2 
range) the dependence on #13 is mild: sin 2 #23 can reach 10~ 6 for 613 ~ 1° or 6 x 10~ 6 
for 0i3 ~ 60°. 

• Sensitivity to sin 2 #13: ji" disappearance 

In Fig. f| we present the 90 % CL exclusion curve in the sin 2 6*13/ Am| 3 plane, at 
different values of #23 = 1°, 10° and 30°, for the \x~ disappearance channel. In |36| it 
was observed that this channel proves more sensitive to sin 2 13 than the /i + appearance 
one for small values of # 23 . On the contrary, the /i + appearance channel has the larger 
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Figure 3: Sensitivity reach in the sin 2 623/ Am^ P^ ane a t different values of 9\% = 
1°, 10°, 30° and 60° for fi + appearance in the 2+2 scheme. 



> 



10" 1 



10" 2 




(>■_■. = 30° 
<r : , = 1°,10° 



10" 
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Figure 4: Sensitivity reach in the sin 2 ^ 13 /Am| 3 plane at different values of6 2 -z 
and 30° for \i~ disappearance in the 2+2 scheme. 



1°,10° 



sensitivity attained for large values of #23; a scenario somewhat disfavoured by the 
LSND measurement. In the /i~ disappearance channel, the Neutrino Factory can put 
an upper bound to sin 2 # 13 at the 10~ 4 — 10~ 2 level for Am^ in the 10 _1 — 10 1 eV 2 
range. 
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Figure 5: Sensitivity reach in the sin 2 ^ 14 /Am| 3 plane at different values of8 13 = 1°, 10° 
and 30° for r~ appearance in the 2+2 scheme. 

• Sensitivity to sin 2 6u and sin 2 #24: r~ appearance 

The t~ appearance channel is quite sensitive to both sin 2 #14 and sin 2 #24- Fig. [| 
illustrates the sensitivity to sin 2 #14 as a function of #13: for about 1°, sensitivities of the 
order of 10 -2 are attainable, while for 10° values as small as 4 x 10 -5 can be reached. 
For even larger values of #13 it goes down to IO -6 (we recall that #13 is not severely 
constrained by the Bugey-CHOOZ experimental bounds, eq. (p6|)). 

Fig. |B| depicts the foreseeable sensitivity reach to sin 2 6*24 as a function of #23 : for 
small values of #23 the sensitivity to sin 2 #24 attains level as low as 10~ 6 . 

In contrast, the r + appearance channel looks less promising, for £3 = 0. Due to the 
relative negative sign between the two terms in the analytic expression for P{y e — * u T ), 
eq. (0), cancellations for particular values of the angles occur, resulting in a decreasing 
sensitivity in specific regions of the parameter space. This sensitivity suppression is 
absent in the r~ channel as the relative sign between the two terms in Piy^ — ► u T ), 
eq. (|8p, is positive^. 

3.2 Sensitivity in the 3+1 scheme 

In the one-mass dominance approximation, the CP-conserving parameter space consists 
of three rotation angles (6m, 624 and #34) and no phases. The useful channels to measure 

9 The same argument holds, albeit interchanging the two channels, for ^3 = n. 
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Figure 6: Sensitivity reach in the sin 2 #24/Am 2 3 P^ ane a t different values of '6 2 z = 1°, 10° 
and 30° for r~ appearance in the 2+2 scheme. 

or put severe upper limits on the gap-crossing angles at the Neutrino Factory are (for 
a ii~ decay): 

v e &n — > f> + (n + appearance) 

V/j, — ¥ v T — > t~ (reappearance). (40) 

We will see in the following that these two channels are optimal to study the whole 
CP-conserving 3+1 parameter space. 

We adopt the same approach as for the 2+2 scheme: we vary sin 2 9 between 10~ 7 
and 1; the other two angles are considered to be already known: one of them is fixed 
to a small value, % = 2°, and the second one is varied from 1° to 60°. The remaining 
parameters (those measured in solar and atmospheric experiments) are taken, following 
H,as: 

9 12 = 22.5°, 013 = 13°, #23 = 45°; 
Am 2 2 = l(T 4 eV 2 , Am 2 3 = 3.5 x l(T 3 eV 2 . 

We consider 2 x 10 20 useful muons per year and 5 years of data taking, with = 20 
GeV . 

• Sensitivity to sin 2 6*14 and sin 2 # 24 : /i + appearance 

The /i + appearance channel is particularly sensitive to both sin 2 # 14 and sin 2 # 2 4- 
Fig. [7| shows the 90 % CL exclusion curve in the sin 2 6 , 14 /Am 2 4 plane (left) and in the 



18 



io- 5 



io- 4 



l(T J ICr 1 
sinVl4) 



IO" 1 



to- 5 



to- 4 



l(T J ICr 1 



IO" 1 



Figure 7: Sensitivity reach in the sin 2 6*14/ Am^ p/ane (7e/£,) and m £/ie sin 2 6*24/ Am 
p/ane (right) at different values of #34 
scheme. 
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1°, 10°, 30° and 60° /or /i + appearance in the 



sin 2 6*24 / Am| 4 plane (right) for different values of #34. The dependence on #34 is very 
mild for small values of 634. In the LSND-allowed region, Am§ 4 6 [10 -1 , 10 1 ] eV 2 , 
both sin 2 14 and sin 2 # 2 4 can reach 10" 4 for 6 34 < 30° or 1(T 3 for # 34 ~ 60°. 

The n~ disappearance channel is not sensitive to #14, but can explore approxi- 
matively the same region as the appearance channel in the sin 2 ^/Am 2 ^ plane, see 
Fig.|. 



to 1 

"a UP 




UT 3 IO" 4 10" 3 LO" 1 to -1 L0° 



sin 2 ) 

Figure 8: Sensitivity reach in the sin 2 ^ 2 4/Am| 4 plane at different values of 634 = 
1°, 10°, 30° and 60° for jjT disappearance in the 3+1 scheme. 
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Figure 9: Sensitivity reach in the sin 2 ^/Am^ plane at different values of 9u = 
1°, 10°, 30° and 60° for r + appearance in the 3+1 scheme (left) and at different values 
°f — 1°, 10°, 30° and 60° for r~ appearance in the 3+1 scheme (right). 



• Sensitivity to sin 2 #34: r appearance 

Both r appearance channel are equally sensitive to sin 2 #34, as can be seen in eqs. 
(|T7|) and (p~8p. Fig. |9| shows the 90 % CL exclusion curve in the sin 2 6*34/ Am| 4 plane 
for t + appearance at different values of #14 (left) and for r~ appearance at different 
values of (9 24 (right). In the LSND-allowed region, sin 2 #34 can reach some units in 10~ 5 
for 0i4, # 24 < 30° or 10~ 5 for 14 , 9 2i ~ 60°. 

3.3 Remarks and conclusions on the sensitivity reach 

The results of the previous subsections show that a Neutrino Factory with n M = 2x 10 20 
useful muons per year and a small detector of 0(1) ton size with r tracking and (/i, r) 
charge identification capability can severely constrain the whole four-family model CP- 
conserving parameter space, both in the 2+2 scheme and 3+1 scheme. In the former, 
the sensitivity reach to all gap-crossing angles in the LSND-allowed region is at the 
level of sin 2 9 > 10 -6 — 10~ 4 , depending on the specific angle considered. In the latter 
the sensitivity reach is at the level of sin 2 9 > 10~ 5 — 10 -3 , depending on the specific 
angle considered, slightly less than in the 2+2 case. 

This results can be easily understood in terms of a simple power counting argument. 
Consider the gap-crossing angles equally small, sin^ ~ e (with 9ij some gap-crossing 
angle in the 2+2 or the 3+1 scheme). In the 2+2 scheme, the CP-conserving transition 
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probabilities become: 



p2+2 
r CP 


> e - 




= 4e 2 sin 2 


p2+2 
r CP 




-> I/ r ) 


= 4e 2 sin 2 


p2+2/ 
r CP \ 




-> z/ T ) 


= 0(e 4 ), 
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-> ^) 
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V / 
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\ 4£ 



+ 0(e 4 



We notice that, with the exception of z/ M 



2A- 



the transition probabilities are generi- 



cally of 0(e ). In the 3+1 scheme, on the contrary, 
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all the appearance transition probabilities are generically of 0(e 4 ). This explains the 
(slight) decrease in the sensitivity in the 3+1 scheme with respect to the 2+2 scheme. 
The 2+2 Up — > v T case is similar to the generic situation in the 3+1 scheme: Fig. ^ 
and Fig. |9] (right) show the same sensitivity reach, indeed. 

Finally, we also present the MNS mixing matrix in the two schemes at 0(e): 



U 
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V -6 
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V-e -6 -e 


1 / 



0(6 2 ) 



(41) 



+ 0(e 2 



(42) 



We remind that in the 2+2 scheme the sterile neutrino is in the first row, v a = 
{i/ s , v e , i/p, i/ T }, whereas in the 3+1 scheme is in the last one, v a = {u e , z/ M , u T , u s }. 
We can build the main contributions to the transition probabilities first writing i/ a , vp 
as linear combinations of mass eigenstates with coefficients given in eqs. (41) and (f42[), 
and then computing \(u a (t)\up)\ 2 . In this way it is simple to derive the behaviour of 
all the transition probabilities P{y a — > vp). 



4 CP-violating Observables 

Genuine CP-violating effects manifest only when at least two mass differences are si- 
multaneously non-vanishing. In the three-family model, the CP-violating contribution 
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to the oscillation probabilities can be written as 



IBS 



Pop = ±2 J ( sin 2Ai2 + sin 2A23 — sin 2A13) (43) 

with J = Ci 2 Si2cf 3 Si3C 2 3S 2 3 sin5 the Jarlskog factor and Ay as defined in eq. (|3]) (the 
± sign refers to neutrinos/antineutrinos). If Ai 2 <C A 23 , Pqp is negligible. Therefore, 
for three-family neutrino mixing the size of the CP-violating oscillation probability 
depends on the range of AmJ 2 , the solar mass difference. In |3B|, |33J it has been shown 
that a maximal phase, \8\ = 90°, can be measured at 90% CL if the LMA-MSW 
solution with Am\ 2 > 2 x 10~ 5 eV 2 is considered. For smaller values of the solar mass 
difference, it seems impossible to measure 5 with the foreseeable beams. However, in the 
four- family model the situation is totally different [|57| : we can consider CP- violating 



observables that do not depend on A so i, but on A atm and Alsnd only. Therefore, 
for four-family neutrino mixing, large CP-violating effects are possible (depending on 
the specific value of the phases 81). In the two-mass dominance approximation the 
parameter space consists of five rotation angles and two phases, both for the 2+2 and 
3+1 schemes. In the 2+2 scheme, the CP-violating oscillation probabilities are given 
by eqs. ( |TT| - p~3D . We notice that, in these expressions, the size of the CP- violating 
probability is linearly dependent on the atmospheric mass difference, A 34 , whereas 
the location of the maximum depends on the LSND mass difference, A 23 . Therefore, 
we expect a maximum in the CP-violating observable at O(10) Km for neutrinos of 
E v = O(10) GeV. With such a short baseline, matter effects are completely negligible. 
In the 3+1 scheme, eqs. ( E?T| - f23| ) , similar results are obtained for A atm = A 23 and 
Alsnd = A 34 . This has to be compared with the three-family model, where the size 
of the CP-violating probability depends linearly on A so ; and the maximum location 
depends on A atm : in this case, the maximum of the CP-observable is expected at 
0(1000) Km, and therefore matter effects are extremely important |58f . 



CP-odd effects are observable in appearance channels, while disappearance ones 
are only sensitive to the CP-even part. The easiest way to measure CP- violation in 
oscillation is to build a CP-asymmetry or a T-asymmetry |]30|] : 



P{v a -> Vp) - P{v a -> v p ) 



a/3 P{v a ^vp)+P{v a ^vp) ' 1 1 



P{v a -> Vp) - P{vp -> v a 



a!i ~ P(u a ^ up) + P(up ^ u a ) ■ 1 ] 

A^p and A^p are theoretically equivalent in vacuum due to CPT, and matter effects 
are negligible at the short distances under consideration. Their extraction from data at 
a Neutrino factory is quite different, though. Consider, as an example, the [y e — > v^) 
channel. The CP-asymmetry, A^, would be measured by first extracting P(u e — > v^) 
from the produced (wrong-sign) yU~s in a beam from fi + decay and P{p e — > v^) from the 
charge conjugate beam and process. Notice that even if the fluxes are very well known, 
this requires a good knowledge of the cross section ratio crip^ — ► /i + )/a(z/ A1 — ► 
Conversely, the measurement of the T-asymmetry, A^ , requires to consider Piy^ — > 
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u e ) and thus a good e charge identification, that seems harder to achieve from the 
experimental point of view. In the following we will deal only with CP-asymmetries. 

A central question on the observability of CP-violation is that of statistics. We do 
not exploit here the explicit E v dependence of the CP-odd effect, and we consider the 
neutrino-energy integrated quantity: 



,cp m = {^^]/iv [/-]}-{iv[U]/iv K]} 

a " [ } {N[l- fi ]/N [l-}} + {N[l+}/N [l+}} ' 



where l a , lp are the charged leptons produced via CC interactions by u a , z/g, respectively 
(the sign of the decaying muons is indicated by an upper index). N[Vg] is the number 
of CC interactions due to oscillated neutrinos, whereas iV [7*] is the expected number 
of CC interactions in the absence of oscillations. In order to quantify the significance 
of the signal, we compare the value of the integrated asymmetry with its error, AA^p, 
in which we include the statistical error and a conservative background estimate at the 
level of 10" 5 . 

In what follows we used the exact expression for the oscillation probabilities, thus 
including the small A so i mass difference and the matter effects. The irrelevance of the 
latter at the considered baseline can be seen in Fig. [IT], where A^ (5 = 90°)/AA^f in 
the 3+1 scheme is presented: we see that matter effects start to be relevant at 0(1000) 
Km. 

Since the matter effects are negligible, the scaling laws with the muon energy 



and the baseline L of the signal-to-noise ratio of the CP- violating asymmetry in eq. (f46|) 
are equivalent to those obtained in vacuum, 
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4.1 CP-violation in the 2+2 scheme 



(47) 



We recall here the results of fl36| , |37|, albeit rederived with slightly different input 
parameters. In the conservative assumption of small gap-crossing angles, #13, #14, #23 
and #24, we consider the following values for the parameters of the MNS mixing matrix 
in the two-mass dominance approximation: 

#13 = #14 = #23 = #24 = 2° ; 0i2 = 45°, #34 = 45° ; 
Am\ 2 = 10~ 4 eV 2 , Am 2 4 = 3.5 x 10~ 3 eV 2 , Am 2 3 = 1 eV 2 ; (48) 

eV 2 



A = 1.1 x 10 



-4 



GeV 



The detector characteristics have been given in Sect. [2T2 . 
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Figure 10: Signal-to-noise ratio of the CP-violating asymmetry in the v e — > channel 
with and without matter effects, for = 10,20 and 50 GeV, as a function of the 
baseline L. The parameters are: Am^ ( = 10~ 4 eV 2 ; Am 2 tm = 3.5 x 10~ 3 eV 2 
/ ^Isnd = 1 eV 2 ; 12 = 22.5°, £ 13 = 13°,0 23 = 45°; 6 U = 6 24 = 6 U = 5°; 
5 2 = o ,#3 = 90°. The matter parameter, A = 2E U \/2Gpn e (with n e the electron 
density in the Earth) is taken to be constant, A — 1.1 x 10~ 4 eV 2 / GeV. This value 
is consistent with a baseline completely contained in the Earth crust ]5S|/ , true for 
L < 4000 Km. 



In Fig. 11 we show the signal-to-noise ratio of the integrated CP asymmetry, eq. 



(flq), in the v e — > channel, as a function of the distance L for three values of the 
parent muon energy, = 10, 20 and 50 GeV. Matter effects, although negligible, are 
included. For this reason, we subtract to the total asymmetry A^f (90°) the matter- 
induced asymmetry, ^4^(0°). We notice that a sizeable signal can be reached: for 
E^ = 50 GeV, approximately 10 standard deviations (sd) at L ~ 30 Km can be 
attained. The scaling of the maximum height with the parent muon energy follows eq. 
(|47l) as expected, increasing with \J~E~ hl . 

• Integrated asymmetry in the v e — > v T channel 



In Fig. [12] we show the signal-to-noise ratio of the subtracted integrated CP asym- 
metry in the v e — > v T channel. The results are pretty similar to those for the v e — > 
channel, with a slightly smaller significance at the maximum: for E^ = 50 GeV , ~ 8 
sd can be attained at L ~ 30 Km. 

• Integrated asymmetry in the v„ — > v T channel 



24 




100 



Figure 11: Signal over statistical uncertainty for CP violation in the v e — > z/ M channel 
in the 2+2 scheme, as a function of the baseline L, for three values of the parent muon 
energy, = 10,20 and 50 GeV. The parameters have been choosen as in eq. ($3j), 
with 5 2 = 0°,5 3 = 90°. 
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Figure 12: Signal over statistical uncertainty for CP violation in the u e — > v T channel 
in the 2+2 scheme, as a function of the baseline L, for three values of the parent muon 



energy, = 10,20 and 50 GeV. The parameters have been choosen as in eq. t\4g), 



with So 



o°,s 3 



90°. 



In Fig. [TBI we show the signal-to-noise ratio of the subtracted integrated CP asym- 
v T channel. The results are totally different from those relative to 



metry in the 



for E,, = 50 GeV, ~ 90 sd can be attained at L ~ 30 Km. 
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Figure 13: Signal over statistical uncertainty for CP violation in the — ► v T channel 
in the 2+2 scheme, as a function of the baseline L, for three values of the parent muon 
energy, = 10,20 and 50 GeV. The parameters have been choosen as in eq. (J&), 
with 5 2 = 0°,5 3 = 90°. 



4.2 CP-violation in the 3+1 scheme 

In the conservative assumption of small gap-crossing angles, #14, #24 and #34, we consider 
the following values for the parameters of the MNS mixing matrix in the two-mass 
dominance approximation: 

#14 = #24 = #34 = 2° ; #12 = 22.5°, #13 = 13°, #23 = 45° ; 
Am 2 2 = 10~ 4 eV 2 , Am 2 3 = 3.5 x 10~ 3 eV 2 , Am 2 4 = 1 eV 2 ; (49) 

eV 2 

A = 1.1 xlO" 4 ^—. 

GeV 

• Integrated asymmetry in the v e — > Vn channel 



In Fig. [L4| we show the signal-to-noise ratio of the subtracted integrated CP asym- 
metry in the v e — > channel, as a function of the distance L for three values of the 
parent muon energy, = 10, 20 and 50 GeV. The results are quite similar to those 
obtained in the 2+2 scheme, although the significance is slightly less: for E^ = 50 
GeV, ~ 6 sd at L ~ 40 Km can be attained. 



The most unfortunate case is presented in Fig. [15]. We show the signal-to-noise ratio 
of the subtracted integrated CP asymmetry as a function of the distance L for three 
values of the CP-violating phases, 62 = S3 = 15°, 45° and 90°. In the 3+1 scheme the 
corresponding oscillation probability, eq. (0), contains three different terms in sin<52, 
sin^3 (with opposite signs) and sin(<52 — #3). Therefore, when 62 = S3 a cancellation 
occurs. This is the reason of the decrease in the significance for ^2 = ^3 = 90° with 
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Figure 14: Signal over statistical uncertainty for CP violation in the u e — > channel 
in the 3+1 scheme, as a function of the baseline L, for three values of the parent muon 
energy, = 10,20 and 50 GeV. The parameters have been choosen as in eq. $4$), 
with 5 2 = 0°,5 3 = 90°. 



respect to the corresponding curve in Fig. [14|. For the lowest value of the phase only 
~ 0.02 sd can be attained. 
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Figure 15: Signal over statistical uncertainty for CP violation in the v e 
in the 3+1 scheme, as a function of the baseline L, for three values of the phases, 
82 = S3 = 15°, 45° and 90°. The parameters have been choosen as in eq. fclQ), with 



= 20 GeV. 
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• Integrated asymmetry in the u e — > u T channel 



In Fig. [16] we show the signal-to-noise ratio of the subtracted integrated CP asym- 
metry in the v e — > v T channel. The results are pretty similar to those for the v e — > 
channel, with a slightly smaller significance at the maximum: for = 50 GeV, ~ 5 
sd can be attained at L ~ 40 Km. 




100 



Figure 16: Signal over statistical uncertainty for CP violation in the u e — > v T channel 
in the 3+1 scheme, as a function of the baseline L, for three values of the parent muon 
energy, = 10,20 and 50 GeV. The parameters have been choosen as in eq. $9Q, 
with 5 2 = 0°,<5 3 = 90°. 



• Integrated asymmetry in the — > v T channel 



In Fig. 17 we show the signal-to-noise ratio of the subtracted integrated CP asym- 
metry in the — > v T channel. Again, as in the 2+2 scheme, the results are totally 
different from those relative to v e — > v^, v T . For E^ = 50 GeV, ~ 100 sd can be attained 
at L ~ 40 Km. 



In Fig. [18| we show the signal-to-noise ratio of the subtracted integrated CP asym- 
metry as a function of the distance L for three values of the CP-violating phases, 
^2 = ^3 = 15°, 45° and 90°. The oscillation probability at the leading order only de- 
pends on sin<5 2 , see eq. (^3]). In this case, for the lowest value of the phase ~ 10 sd 
can still be reached. 



4.3 Remarks and conclusions on the CP-violation observables 

The previous results clearly show how a maximal CP-violation is easily measurable 
with a not-so-large detector of 10 Kton size, at a baseline L = O(10) Km. The optimal 
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Figure 17: Signal over statistical uncertainty for CP violation in the — ► v T channel 
in the 3+1 scheme, as a function of the baseline L, for three values of the parent muon 
energy, = 10,20 and 50 GeV. The parameters have been choosen as in eq. (J&), 
with 5 2 = 90°, 5 3 = 0°. 




Figure 18: Signal over statistical uncertainty for CP violation in the — > v T channel 
in the 3+1 scheme, as a function of the baseline L, for three values of the phases, 
^2 = ^3 = 15°, 45° and 90°. The parameters have been choosen as in eq. with 



= 20 GeV. 



channel to observe CP- violation is the — > v T channel, where 0(100) sd can be 
attained. Even at lower energies, ~ 40 — 50 sd can be easily found. The number 
of expected A^ T ± in the detector is of O(10 5 ). Notice that in the z/ M — > v T channel 
a non-maximal CP-violating phase gives a significance at the level of ~ 10 sd even 
for the lowest energy, E^ = 10 GeV. The other two channels, u e z/ M , v T are quite 
similar and give much smaller significance. For a smaller detector mass M, a reduction 
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factor oc 1/ y M should be applied to all the results of this section. Therefore, for an 
OPERA-like 1 Kton detector we still expect large CP-violating effects in the — ► v T 
channel. 

All of these results are obtained in both the 2+2 and the 3+1 scheme, with slight 
differences between the two. The real gain with respect to the three-family model is 
that the small solar mass difference, that modules the overall size of the CP-violating 
asymmetry, is traded with the much larger atmospheric mass difference. 

In the three-family model, the largest CP asymmetry is expected in the v e — ► 
transition. We have shown that in the four-family model (both in the 2+2 and 3+1 
scheme) this is not the case: it is — > v T the optimal channel to study CP violation. 
The reason can be easily understood applying the power counting argument introduced 
in Sect. p75L Neglecting the background, the signal-to-noise ratio of the CP asymmetry 
AST/AA^is: 

A °^ cp K p ar(y« -> u p) (50) 

aa *p \jP C p{y<x -»• vp) 

The CP-violating transition probabilities in the approximation of equally small gap- 
crossing angles, sin 6^- ~ e (with 9%j a generic gap-crossing angle in the 2+2 or the 3+1 
scheme), in the 2+2 scheme become: 



^Ve - v») = "4e 2 sin(2fl 34 ) sin(5 2 + 6 3 ) ( ] sin 2 ( ] + 0(e^ 



( Amj 4 L \ 2 / Am 2 3 l / 
I AE V ) I AE V 



(51) 



- W) = +^ M^u) sin(5 2 + (^fy sin 2 (^^) + 0(e 4 ) 



(52) 



P&\vh^Vt) = -4e 2 sin(2M[sm5 2 + sin(5 2 + 5 3 )] f^^) sin 2 (^^^+0(e 4 ) 



(53) 



in the 3+1 scheme become: 

F&ly^v*) = 4 e 3 sin(2Msin(5 2 -5 3 ) ^(^^]+0(e^ 



V 



^V- - *r) = -8e 3 c 2 3 sin5 3 (^J^) sin 2 + 0(e 4 ) , 



^%^) = -4 e 2 sin(2Msin5 2 (^^) sin 2 f^^)+0( e 4 ) 



(54) 



(55) 



(56) 
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Scheme 


Transition 


Pcp 




A/ A A 




v e — > 


e 2 


e 


0(1) 


Three-family 


v e — > V T 


e 2 


e 


0(1) 




VfM -> 


1 


e 


0(6) 




-> ^ M 


e 2 


e 2 


0(e) 


2+2 


Z/ e -> I/ T 


e 2 


e 2 


0(e) 




^ -> 


e 4 


e 2 


0(1) 




-> ^ M 


e 4 


e 3 


0(e) 


3+1 


Z/ e -> I/ T 


e 4 


e 3 


0(e) 




^ -> *+ 


e 4 


e 2 


0(1) 



Table 3: Small angles suppression in the CP-conserving and CP-violating oscillation 
probabilities, and in the signal-to-noise ratio of the CP asymmetries, in the three-family 
model and in both four-family model mass schemes. 



(in this case we have considered S13 ~ e also, since the present experimental results in 
the three-family model show that #13 < 13° ||55|| ). 

In Tab. § we report the leading order in e for the different Pep and Pqp in the 
three-family model and in both schemes of the four-family model. In three-family, the 
small parameter is S13 ~ e. In four-family we consider equally small LSND gap-crossing 
angles: s i3 = s i4 = s 2 3 = s 24 ~ e for the 2+2 scheme; S14 = s 24 = s 34 = e for the 3+1 
scheme. In this last case we take S13 ~ e, also. We also report the leading order in e of 
the signal-to-noise ratio of the various CP asymmetries, using eq. (pd|). 

The last column can be easily read: in the three-family model, the v e — > and 
v e — > v T have a signal-to-noise ratio of the corresponding CP asymmetry of 0(1) in 
the small angles. On the contrary, in both the 2+2 and 3+1 four-family model, it is 
the v T channel to be of 0(1) in the small angles, thus justifying a posteriori our 

results. As a final remark, notice that the last column for the 2+2 and the 3+1 scheme 
is identical, although the corresponding CP-conserving and CP-violating probabilities 
are of different order in the small angles suppression. 



5 A magnetized iron detector with no r-tracking 

In this section we explore the possibility of reconstructing two angles or one angle and 
a phase at a time in a short baseline experiment, L = 40 Km (the distance where 
the CP- violating observable of the previous section is maximized). We focus on the 
—* channel and consider a realistic 10 Kton magnetized iron detector with \x charge 
identification of the type discussed in |8J. A detailed analysis of the backgrounds and 
detection efficiencies of this apparatus can be found in the literature. The Neutrino 



31 



Factory is run with 2 x 10 20 useful muons per year for 5 operational years for both 
muon polarities at = 50 GeV , with a detector energy resolution of AE U = 10 GeV 
. Our motivation is the following: is it possible to avoid the complications connected 
with the r detection taking full advantage of the energy dependence in the v e — ► 
channel, instead? 



We follow the procedure described in |33|: let N ip be the total number of wrong- 
sign muons detected when the Neutrino Factory is run in polarity p = p + , pT , grouped 
in 5 energy bins specified by i — 1 to 5. In order to simulate a typical experimental 
situation we generate a set of "data" Ui tP as follows: for a given value of the oscillation 
parameters, the expected number of events, N i)P , is computed; taking into account 
backgrounds and detection efficiencies per bin, b ijP and e i;P (as quoted in f3~3fl), we 
perform a Gaussian (or Poisson, depending on the number of events) smearing to 
mimic the statistical uncertainty: 

Smear(iV ij ,e ij p + b i;P ) - b Lp 
n i)P = ^ . (57) 

The "data" are then fitted to the theoretical expectation as a function of the mixing 
matrix parameters under study, using a x 2 minimization: 

* 2 = ??P^) 2 < (58) 

where 6rii iP is the error on n^ v (we include no error in the efficiencies). 

First consider the simultaneous measurement of # 14 and #24- In the one-mass domi- 
nance approximation the v e — > transition probability, eq. (|IED, depends on the com- 
bination s 2 4 S2 4 . This term (symmetric under 6*14 <-> #24) dominates over the sub-leading 
A m 2 atm - dependent (non-symmetric) corrections. Therefore, the energy dependence of 
N itP is not enough to resolve the two angles. 

Next, we consider the simultaneous measurement of #34 and one of 814,624,. We 
compute the leading corrections to eq. ([!(]) in powers of Am 2 atm . For vanishing phases 



and S13 = s 2 4 = e , #23 = 45°, the leading terms are 

P^pVe -> Up) = 4e 2 s 2 4 C34 sin 2 A 34 - 2C34S34S14 e A 23 sin(2A 34 ) 

+ 0( e 4 )+0(A 2 3 )+0( e 2 A23). (59) 

For vanishing phases and S13 — su = e, 9 2 3 = 45°, the leading terms are, instead, 

Pcpiy* ^ Vy) = 4e 2 c 2 4S 2 4 4 sin 2 A 34 + 0(e 4 ) + 0(A 2 3 ) + 0(e 2 A 23 ) • (60) 

In Fig. |19| (left) we present the 68, 90 and 99 % confidence level contours for a si- 
multaneous fit of 9 u and #34 for a "data" set generated with 9 14 = 5°, #24 = 2° and 
#34 = 6°. This figure corresponds to eq. (|59|), where the leading correction to the one- 
mass dominance formula is 0(eA 2 3). In Fig. |l^ (right) we present the 68, 90 and 99 % 
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Figure 19: 68, 90 and 99 % CL contours resulting from a \ 2 fit of 9 U and #34 (left) or 
9 2 4 and #34 (right). The parameters used to generate the "data" are depicted by a star: 
6*i4 = 5°, 6*34 = 6° (left); # 24 = 5°, 6*34 = 6° (right). Only statistical errors are included. 



confidence level contours for a simultaneous fit of 624 and #34 for a "data" set generated 
with #24 = 5°, B44 = 2° and # 34 = 6°. This figure corresponds to eq. flS0"|): notice that 
the 0(eA 23 ) correction to the one-mass dominance formula is absent and the leading 
corrections start at higher orders. The sensitivity to #34 is therefore suppressed with 
respect to eq. ((59"|). In summary, # 14 or 9 2 4 are reconstructed with a precision of tenths 
of degree; on the contrary, # 34 is measurable with a very poor precision in both cases, 
with slightly better results in Fig. |19| (left). 

We try also to simultaneously fit the combination S14S24 and #34. In Fig. |(] we 
generate the "data" for Si 4 s 2 4 = 0.09, # 34 = 5°. As in Fig. [T5], we observe that #34 
is poorly reconstructed, whereas Si 4 s 2 4 is severely constrained. Our results seem to 
indicate that is very difficult to measure #34 using the energy dependence of sub-leading 
effects in the v e — > channel. 

Finally, we consider the simultaneous measurement of one gap-crossing angles (the 
other two being fixed at some small value) and one CP-violating phase. In this case, 
the leading correction to the one- mass dominance formula, eq. fll6|) , is the CP-odd 
contribution, eq. (f2~lj). In Fig. |21] we show the confidence level contours for a simul- 
taneous fit of #14 and £3. The theoretical values for which the "experimental data" 
have been generated are: #14 = 3°,5 3 = 50° (6* 2 4 = #34 = 2°). We see that the angle 
is reconstructed with a precision of tenths of degree, whereas the phase is measured 
with a precision of tens of degrees only. This is precisely the same situation of the 
three-family results of f3~3fl , however. 
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Figure 20: 68, 90 and 99 % CL contours resulting from a y 2 fit of S14S24 and #34. The 
parameters used to generate the "data" are depicted by a star, Si 4 s 24 = 0.09, #34 = 5°. 
Only statistical errors are included. 
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Figure 21: 68, 90 and 99 % CL contours resulting from a \ 2 fit of 9 U and 63 (left). 
The parameters used to generate the "data" are depicted by a star, 9u = 3°, 63 = 50°. 
Only statistical errors are included. 



6 Conclusions 

The ensemble of solar, atmospheric and LSND neutrino data can be explained with 
three active plus one sterile flavour states. Although the four-family neutrino mass 
spectrum preferred by the experimental data is the so-called 2+2 scheme (with two 
almost degenerate pairs well separated by the LSND mass difference), the latest LSND 
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result are marginally (at 99 % CL) compatible with the 3+1 scheme (three almost 
degenerate neutrinos separated from a fourth, mainly sterile, one). 

In this paper we studied the physical potential of a Neutrino Factory in both four- 
family model schemes, in the spirit of |36| , |37[| where the 2+2 scheme was carefully 
examined. A re-analysis of the 2+2 results and a totally novel analysis in the 3+1 
scheme has been presented. We have derived one- and two-mass scale dominance 
approximations appropriate for CP-even and CP-odd observables, respectively, in both 
schemes. 

The rich flavour content of a muon-decay based beam is extremely useful to deter- 
mine or severely constrain the four-family model parameter space: in both schemes, 
the sensitivity to gap-crossing angles as small as sin 2 6^ = 1CT 6 — 1CT 4 (depending on 
the specific angle) can be achieved, with a 1 ton detector at L ~ 1 Km down from the 
source, for n (i = 2 x 10 20 useful muons per year and 5 years of data taking. In the 3+1 
scheme, we notice a slight loss in sensitivity with respect to the 2+2 scheme, that we 
interpret as a consequence of the higher power in the small gap-crossing angles in the 
leading terms of the CP-conserving transition probabilities. 

CP violation may be easily at reach with a 10 Kton detector at L = O(10) Km, 
especially through "r appearance" signals, in both schemes. The increased significance 
of the CP-violating observables with respect to the three-family model asymmetries 
is due to the fact that asymmetries proportional to Am 2 tm are possible, whereas in 
the three-family model the CP-violating asymmetries are proportional to the small 
Am^j. Moreover, in the four-family case the asymmetries are modulated by Am'j jSND , 
thus peaking at L — O(10) Km; in three families they are modulated by Am 2 tm and 
therefore they peak at L — 0(1000) Km, thus significantly suffering from matter effects. 
On the contrary, matter effects are totally negligible in the four-family case. We give a 
simple argument based again on the power counting of the small gap-crossing mixing 
angle to explain why in four-family models (both in the 2+2 and 3+1 schemes) the 
v p, ~^ v t channel is the optimal one to study CP- violating observables, whereas in three 
families u e — > seems to be best suited. 

Eventually, we consider a 10 Kton detector at L = 40 Km with only \i charge iden- 
tification (of the magnetized iron type). In this case, we use the energy dependence of 
the wrong-sign muons to try to simultaneously reconstruct two angles (or combinations 
of them) or a CP-violating phase and one angle at a time. Our results suggest that it is 
not possible to constrain the whole CP-conserving parameter space using the v e — > 
channel only. However, we can simultaneously measure one angle and one phase, with 
a precision of tens of degrees in the latter. 

Summarizing, a Neutrino Factory has an enormous discovery potential when four- 
neutrino models are considered, in both the 2+2 and 3+1 schemes. If the LSND results 
will be confirmed by MiniBooNE, a muon storage ring appears to be an extremely 
powerful facility to perform precision measurement on the whole four-family model 
parameter space and most probably the best opportunity to discover CP violation 
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in the leptonic sector. A (not-so-large) detector with r-tracking and (//, r) charge 
identification capability is needed. If MiniBooNE will not confirm LSND, however, 
the 3+1 scheme will still represent a possible extension of the Standard Model and a 
Neutrino Factory can severely constrain its parameter space. This cannot be said of 
the 2+2 scheme, that would be ruled out by a negative result of MiniBooNE. 
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